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In this article we obtain a fundamental measure functional for the model of aligned hard hexagons in the 
plane. Our aim is not just to provide a functional for a new, admittedly academic, model, but to investigate the 
structure of fundamental measure theory. A model of aligned hard hexagons has similarities with the hard disk 
model. Both share "lost cases", i.e. admit configurations of three particles in which there is pairwise overlap 
but not triple overlap. These configurations are known to be problematic for fundamental measure functionals, 
which are not able to capture their contribution correctly. This failure lies in the inability of these functionals 
to yield a correct low density limit of the third order direct correlation function. Here we derive the functional 
by projecting aligned hard cubes on the plane x + y + z — 0. The correct dimensional crossover behavior of 
these functionals permits us to follow this strategy. The functional of aligned hard cubes, however, does not 
have lost cases, so neither had the resulting functional for aligned hard hexagons. The latter exhibits, in fact, 
a peculiar structure as compared to the one for hard disks. It depends on a uniparametric family of weighted 
densities through a new term not appearing in the functional for hard disks. Apart from studying the freezing 
of this system, we discuss the implications of the functional structure for new developments of fundamental 
measure theory. 
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I. INTRODUCTION 

Fundamental measure (FM) theory [ll] is one of the 
most successful density functional (DF) theories, yet 
also one of the most difficult to adapt to new sys- 
tems. Unlike classical DF approximations |2], which 
describe a general approximate recipe in which some 
knowledge of the fluid is "cooked up" to produce 
a functional, FM theory constructs functionals from 
geometric principles in a far more involved manner. 
As a result of this process the resulting functionals 
have got nicer features; among them, three are strik- 
ing: their higher predictive power (they yield structure 
functions that are needed as input in classical approx- 
imations), their natural formulation for multicompo- 
nent mixtures, and their good behavior under dimen- 
sional crossover (when d-dimensional systems are con- 
strained to d — 1 dimensions, rf-dimensional FM func- 
tionals become (d — 1) -dimensional ones). However, 
the price to pay for having such nice functionals is that 
their structure is extremely rigid: almost any reason- 
able modification one makes to "improve" the quality 
of the results spoils one of the above features, mainly 
dimensional crossover JH 0, HI. The latter is not only 
a remarkable property that FM functionals (and only 
them) share with exact ones, but also a desirable prop- 
erty of any functional which is meant to study fluids 
under strong confinement IH 0,111]. 

This is the reason why every extension of FM theory 
beyond the hard spheres fluid for which it was origi- 
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nally proposed ifHH [Toll has become a "major achieve- 
ment". Extensions are nowadays available for parallel 
hard cubes and parallelepipeds ifTTl [l2l [TH [Till (which 
prov ide a restricted orientation model of liquid crystals 
lll5ll ): soft spherical potentials iflfjl : non-additive mix- 
tures J jjl Ta, 19|l: mixtures of rods, plates and spheres 
Jlp-JlIUp El HI; lattice fluids IH J2I 
u9l l30l Bill ; and fluids in porous media 1321 l33l 1341 . 
Even for hard spheres, Rosenfeld's original functional 
|Q]] has undergone important improvements over the 
years ll3^l3^l37ll38ll . after realizing that dimensional 
crossover was a unique feature of this type of function- 
als very much entangled to its construction procedure. 

When Rosenfeld first conceived FM theory (TJ] it 
rested strongly on geometrical properties of spherical 
overlaps and on scaled-particle theory 13911 . A decade 
later, the theory had been reformulated in terms of 
"zero-dimensional (0D) cavities" IB~5l l36ll . By this it 
must be understood a cavity able to hold no more than 
a hard sphere. Under the requirement that confinement 
of the FM functional to one such cavity must lead to 
the exact result, and introducing one- two- and three- 
point cavities, by adding and subtracting the necessary 
terms so as to maintain the exact 0D limit a functional 
arises with the required structure IB~7l l38ll . The result 
is not "perfect" in the sense that there are three-point 
0D cavities for which the exact result cannot be recov- 
ered: those for which three spheres placed at the three 
points of the cavity have pairwise overlap but not triple 
overlap. These cavities were termed "lost cases" l37ll 
because they do not contribute to the free energy, and 
their existence reveals the inability of this construction 
to reproduce the lowest order in the density expansion 
of the three-particle direct correlation function (DCF) 
Jll] (which is non-zero for those configurations). As 
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a matter of fact, the problems arising from the exten- 
sion of this functional to mixtures of hard spheres have 
the same origin, and corrections trying to palliate these 
problems are unable to remedy the defect of the corre- 
lations 051]. 

On the contrary, the FM functional for aligned hard 
parallelepipeds does not have lost cases because for 
this kind of particles, whenever three particles have 
pairwise overlap there is necessarily triple overlap. In 
fact, the exact OD limit is recovered for cavities of 
any shape [12], and the low-density limit of the three- 
particle DCF is exact. This led to the belief that FM 
theory is simply unable to produce a functional for hard 
spheres without lost cases 

The extension of FM theory to lattice fluids is based 
on this OD cavity reformulation H3 E3, |H |H 
l30l HU]. But because lattice geometry lacks spherical 
symmetry, this extension has uncovered an important 
ingredient in the theory. While the FM functional for 
hard spheres is expressed in terms of a set of weighted 
densities whose weights are associated to geometrical 
features of the particles B E3 IH HI 111 the 
weights in lattice FM functionals are associated to ge- 
ometrical features of maximal OD cavities IT261, 129I l30ll . 
These are OD cavities such that if they get extended 
in any way, they are not OD cavities anymore 114 ill . 
Maximal OD cavities need not have the same shape as 
the particles that define them (in most cases they will 
not); for instance, maximal OD cavities of a nearest- 
neighbor exclusion lattice gas in a triangular lattice, 
which is represented by hard hexagons, are equilat- 
eral triangles (see Ref. [28] for this and other exam- 
ples). At the same time, by construction these lattice 
FM functionals recover the exact OD limit for any OD 
cavity ||29[ l30ll and as a consequence of this fact, they 
can be proven to yield the correct low-density limit of 
the three-particle DCF Il30fl . Interestingly, maximal OD 
cavities for aligned hard parallelepipeds have exactly 
the same shape as the particles, however this is not true 
for spheres, where apart from spherical cavities there 
are other maximal OD cavities with different shapes 
(see Fig.[T]for an example). 

In Ref. [28] FM functionals for many two- and three- 
dimensional hard core lattice gases were obtained from 
the known functional of hard (hyper)cubes in a (hy- 
per)cubic lattice by exploiting dimensional crossover. 
One of these was the hard hexagons model, which was 
obtained from the hard cubes model by constraining 
the centers of mass of the cubes to lay on the plane 
x + y + z = 0. In this paper we will apply the same 
procedure to the continuum FM functional of paral- 
lel hard cubes in order to obtain the (continuum) two- 
dimensional fluid of aligned hard hexagons. The rea- 
son to do this is the following. Three aligned hard 
hexagons can be arranged in such a way that there 
is pairwise overlap but no triple overlap (see Fig. [2]); 
hence, according to the cavity construction of the FM 
functional for hard spheres IB71 l3al . there should be 
"lost cases". However, we are going to obtain such 




FIG. 1: (color online) Example of a non-spherical maximal 
OD cavity for the system of two-dimensional hard spheres 
(disks) of radius R (colored circle in the figure). Any rota- 
tion of this cavity will produce a new maximal OD cavity. To 
illustrate that this cavity is not contained in a spherical maxi- 
mal OD cavity, the latter is plotted on top with dashed line. 




FIG. 2: (color online) Configuration of three aligned 
hexagons exhibiting pairwise, but not triple overlap (lost 
cases). 



a functional by dimensional crossover from the func- 
tional of parallel hard cubes which does not have lost 
cases. As explained, this means that the OD limit is re- 
covered for any OD cavity; therefore the same will hold 
for the resulting functional for aligned hard hexagons. 

The FM functional for aligned hard hexagons that 
we will obtain provides interesting insights into FM 
theory. First of all, the maximal OD cavities, not the 
particles, are the relevant constructive geometrical ob- 
ject. Secondly, it points out that a FM functional for 
hard spheres (or disks) will probably have an infinite 
number of terms. From a practical point of view this is 
good and bad news: good, because we know what the 
FM functional for hard spheres IT38ll is missing in or- 
der to get rid of the lost cases; bad because a functional 
with an infinite number of terms will be useless for real 
purposes. At the end of this paper we will discuss these 
issues in more depth. We think, however, that despite 
its eventual utility, the information that this FM func- 
tional for aligned hard hexagons provides is relevant 
for a thorough understanding of FM theory. 

The rest of the paper is organized as follows. In 
section [TT] we describe the construction of the FM 
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FIG. 3: (color online) Constraining the centers of mass of 
parallel hard cubes to lay on the plane x + y + z — yields 
a system of aligned hard hexagons. 



functional for aligned hard hexagons by dimensional 
crossover of the FM functional for the fluid of paral- 
lel hard cubes. The procedure as well as the resulting 
weighted densities and the form of the functional are 
explained in this section, but the detailed calculations 
are deferred to Appendix [A] In Sec. [Ill] we analyze 
the equation of state derived from the functional both 
for the uniform fluid and for the triangular solid phase. 
The functional predicts a first order freezing at coex- 
isting packing fractions r\f — 0.58 for the fluid and 
T] s = 0.63 for the solid. Finally, in section HVl we dis- 
cuss the features of the resulting functional, with spe- 
cial emphasis on those that cause the functional to be 
free from lost cases. 



II. FROM HARD CUBES TO HARD HEXAGONS 



As described in Ref. 12811 the way to obtain an effec- 
tive system of aligned hard hexagons is to start off from 
a system of hard cubes aligned parallel to the coordi- 
nate axes and constraint their centers of mass to lay on 
the plane x + y + z = 0. Figure [3] illustrates this geo- 
metrical construction. Making use of the good behav- 
ior of FM functionals under any dimensional crossover, 
we will carry out the same projection in the FM func- 
tional of parallel hard cubes and thus obtain the one for 
aligned hexagons. 

The FM functional of parallel hard cubes of edge- 
length 2L can be written as lfl2l[l3"ll 
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j3 = (kT) , with k the Boltzmann constant and T the 



temperature) for a density profile of the hard cube fluid 
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is the excess, over the ideal, free energy. The functions 
p a (r) are weighted densities 



p a (r) = J dr' uj a (r - r')p(r'), 
where the scalar or vectorial weights are given by 



w 3 (r) = T(x)T(y)T(z), 
w 2 (r) = (C{x)r(y)T(z),T{xX(y)r(z), 
T(x)T(y)((z)) , 

«i(r) = (T(a;)C(y)C(*),C(3Mv)C(*), 
C(x)C(y)T(z)), 

wo(r) = C(s)C(l/)Cto, 
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with 

r(«) = e(L-H), 



((u) = 6(L-\u\), (13) 



Q(x) being Heaviside's step function (0 if x < and 
1 if x > 0) and 6(x) Dirac's delta. Notice that p2j 
(J — 1, 2, 3) denotes the jth component of P2. 
Now, the projection amounts to taking 



p(r) =p(x)5(x + y + z) 



(14) 



in the functional (fl}, where p(x) = p(x, y) is 
the density profile of aligned hard hexagons. The 
choice of coordinates corresponds to a change to 
the (non-orthogonal) basis {t>i,b2,b3} given by 
(bi,b 2 ,b 3 ) = (ei,e 2 ,e 3 )P, with 



P 




(15) 



and {ei, e 2 , e 3 } the canonical basis. Vectors bi andb 2 
form a basis on the plane x+y+z = (see Fig.|4j». This 
choice of vectors amounts to working with the projec- 
tions of the hexagons on the XY plane, because the 
projections of bi and b 2 are simply ei and e 2 . 

The details of introducing the density profile ( fT4b 
into the excess part of the free-energy functional are 
deferred to Appendix [A] Here we simply give the final 
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u=L 



0<u<L 




-L<u<0 u=-L 

FIG. 5: (color online) Maximal OD cavities for the system 
of aligned hard hexagons are obtained as the sections of the 
cubes by planes x + y + z + u = with — L < u < L. 
Triangular cavities correspond to u = ±L while hexagonal 
ones to — L < u < L (u = is the regular hexagon). 



result. The projection transforms the original weighted 
densities for the cubes, p a (r), into a set of new densi- 
ties for the hexagons. The most striking result is that 
these new weighted density are associated to maximal 
OD cavities for the hexagons, not to the hexagons them- 
selves. The complete set of such maximal OD cavities 
can be obtained as the sections of one of the origi- 
nal cubes by the planes x + y + z + u = 0, where 
—L < u < L (see Fig. [5]). The cases u = ±L cor- 
respond to two equilateral triangles (pointing up and 
down), while the cases —L<u<L correspond to 
hexagons (of which only u = is a regular hexagon 
identical to the fluid particles). This comes as an im- 
portant difference with respect to the FM functional for 
hard spheres, and in retrospect yields a new interpreta- 
tion to the weighted densities of the original cubes as 
associated to maximal OD cavities (which in the case 
of cubes are indistinguishable from the particles). 

To be precise, the weighted densities that we need to 



FIG. 6: Weighted densities are associated to geometric ele- 
ments of the cavities: rig (x) (a) and rig ' (x) (b) to the ver- 
tices of the triangular cavities; the components of n[ + ' (x, u) 

(c) and n[ ' (x, u) (d) to the edges of the cavities. In (c) and 

(d) the numbers 1 , 2, 3 label the component which is associ- 
ated to each edge. 



express the functional are 

4 ±} (x) = /dx'O^Cx-x'Mx'), (16) 

n[ ±] (x,u) = /dx'nPCx-x'^Mx'), (17) 

na(x,u) = y*dx'r! 2 (x-x',u)p(x'), (18) 

and it is convenient to define also the weighted densi- 
ties 



nP(x) - ni IJ (x,±L) 



(±)/ 



r4 ±J (x) = n, 2 (x, ±L). 



(19) 
(20) 



The weights that define these densities are the follow- 
ing: 

^ ±} (x) = S(xTL)5(y T L) 

+5(x =F L)8(y ± L) (21) 
+S(x±L)8(yTL), 

nf^x.u) = (i(iTl)T(y)T(ti-i-i/), 

T(x)S(yTL)T(u-x-y), (22) 

T(x)r(y)5(u - x - y T L)J , 

f2 2 (x,it) = t{x)t{u)t{u - x - y). (23) 

The meaning of these weighted densities is related to 
averages over different geometric elements of the max- 
imal OD cavities to which they are associated. Thus, 
ri2 (x, u) is the average over the area of the cavity cor- 
responding to that value of u (the colored regions in 
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Fig. |5j; each component of (x, u) is the average 
over one edge of the hexagonal cavity (triangular if 
u = ±L); and (x) is the average over the three 
vertices of the corresponding triangular cavity. The two 
latter cases are illustrated in Fig. [6] 

With the help of these weighted densities we can 
write /3Fjgj H [p] = / dx $ A hh(x), where 
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(For the sake of notational simplicity we have omitted 
the argument x in all weighted densities, retaining only 
the argument u in those that depend on it.) 

There are several features worth noticing in this FM 
functional for aligned hard hexagons which we have 
derived from the one for parallel hard cubes. First of 
all, the most obvious fact: weighted densities are as- 
sociated to the geometry of maximal OD cavities, as in 
lattice FM functionals, and not to the geometry of parti- 
cles, as in the FM functional for hard spheres or disks. 
Secondly, as in the system of aligned hard hexagons 
there is an infinity of maximal OD cavities, the third 
term exhibits a "sum" over them all; hence the integral 
in that term. Finally, the typical FM structure for $ as a 
sum of D terms, D being the dimension of the problem, 
breaks down here: we have a fwo-dimensional system 
which is described as a sum of three terms. 

All these features will have consequences for the 
general structure of FM functionals, which we shall 
discuss later in Sec. [IV] 



given that the area of a hexagon is Vh = i\fiL 2 (re- 
call that L is half the edge length of the cubes), hence 
the packing fraction rj = VhP, the weighted densities 
reduce in this limit to 
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where the packing fraction r\ = v^p, and therefore the 
excess free energy per unit volume (in kT units), $, 
becomes 
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From this expression we can obtain the equation of 
state as 

Pp g(g/p) 

p o-q 



l-r?/2 

(1-77)2 
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Perhaps the most striking feature of this equation of 
state is its divergence as (1 — rj)~ 5 / 2 at close packing. 
This exponent 2.5 is noticeably higher than the expo- 
nent 2 that a straightforward scaled particle argument 
would predict. 



B. Solid phase 



III. THERMODYNAMICS OF THE FLUID OF 
ALIGNED HARD HEXAGONS 

A. Fluid phase 

The free energy of the fluid phase is obtained by spe- 
cializing the weighted densities with a uniform parti- 
cle density. The only subtle point we have to take into 
account is that, because of our choice of coordinates 
(which actually describe the projections of hexagons 
on the XY plane) if p denotes the particle density of 
hexagons (measured on the plane x + y + z = 0), the 
uniform density profile will reduce to p(x) = V3p. 
(This \/3 is the scale factor difference between actual 
hexagons and their projections.) With this in mind, and 



The standard way to approach the solid phase in DF 
theory is to use a parametrization of the density profile 
as a sum of Gaussians centered at the nodes of the solid 
lattice, 

/3(x') = — exp { — a(x' — rda', — sda! 2 ) 2 \. 

ri,r 2 £Z 

(35) 

Here p(x') denotes the density profile of hexagons on 
the plane x + y + z = 0; x' = (x',y') denotes the po- 
sition referred to an orthogonal coordinate system on 
that plane; a is related to the mean square displace- 
ment of the particle with respect to its lattice node; 9 is 
the occupancy of the solid (the mean number of parti- 
cles in one unit cell), which accounts for the vacancies; 
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FIG. 7: Reduced pressure, f3pvh, with Vh = 3\/3Z/ 2 the area 
of a hexagon, versus pacing fraction r\ = pv^ for a fluid of 
aligned hard hexagons. The packing fractions at coexistence 
are rjf = 0.58 for the fluid and rj s = 0.63 for the solid. 



and d = d c ^9/r] is the lattice parameter, with r\ the 
packing fraction and d c the lattice parameter at close 
packing, d c = \[§L (hexagons have edge length V2L). 
The unit vectors ai = (0, 1) and a! 2 = (-• \/3/2, 1/2) 
are a convenient choice for the basis of the triangular 
lattice's unit cell. 

Although the implementation of this density profile 
is rather straightforward, a few words on the appropri- 
ate choice of variables may be helpful. First of all, as 
all weighted densities are expressed in terms of pro- 
jected coordinates on the plane XY, we should describe 
the density profile in terms of these coordinates. Thus 
we can write x' = Jx, with 




(36) 



and, given that dct J = \/3, the density profile be- 
comes p(x) = v / 3p(x / ). In the projected represen- 
tation, the unit cell is a rhombus, so it is convenient 
to introduce in the integrals the change of variables 
x — ► Qx, where 



This transforms the global integral as 



dx 

unit cell 



d/2 rd/2 

dx I dy 

-d/2 J-d/2 



(37) 



(38) 



(notice that dct Q = 1/2). 

In the projected coordinates, the sum of Gaussians 
defining the density profile can be factorized as 



p(x) = \/30 9a +y r2 ) 9a( 
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FIG. 8: Fraction of vacancies, 1 — 8, as a function of the 
packing fraction r\ for the fluid phase. 
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FIG. 9: Mean square displacement of particles around the 
equilibrium positions in the solid phase, as a function of the 
packing fraction 77. 



with the definitions x ri — x — r\d, y r2 = y — r 2 d and 
9a(x) — \Ja/-Ke~ ax . This permits one to express the 
weighted densities as products of Gaussian and error 
functions. 

Minimization of the functional is carried out nu- 
merically in the parameters a and 9. This determines 
the free energy of the solid phase for every packing 
fraction 77. At the point where this free-energy curve 
branches off the one of the fluid phase, the slope dis- 
continuously decreases. This produces a concavity in 
the whole free-energy curve and therefore the transition 
is first order and the coexisting densities are determined 
via a standard Maxwell's double-tangent construction. 
The resulting equation of state, depicting both the fluid 
and the solid pressures, is plotted in Fig. [7] Also plot- 
ted in Figs.[8]and|9]are the fraction of vacancies, 1 — 9, 
and the square root of the mean square displacement of 
hexagons with respect to their lattice positions. Inter- 
estingly, the solid of hard hexagons has a low fraction 
of vacancies all the way up from the transition (never 
larger than 5%), in marked contrast with what happens 
for the fluid of hard cubes from which the functional 
for this system is derived [ 14]. 
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One last remark is in order. The result found here 
for the equation of state of the system of aligned hard 
hexagons is very similar to that of its lattice counterpart 
ll28ll . For the latter the exact result is known to have a 
continuous transition 14011 . although so smooth that a 
first order discontinuity is not a bad quantitative ap- 
proximation. The exact result for the continuum model 
is unknown, but certainly the same caveat on the nature 
of the transition applies to it. 



IV. DISCUSSION 

The fluid of aligned hard hexagons has strong simi- 
larities with the fluid of hard disks. The most important 
for the aims of this work is that particles may be placed 
in configurations such that they have pairwise overlap 
but not triple overlap (see Fig. |2j. These configurations 
have been termed "lost cases" [37] because FM theory, 
as currently formulated for hard spheres or hard disks 
ll37l l38ll . is unable to capture their contribution. The 
basic constructive principle of FM theory is the recov- 
ery of the exact OD limit of the free energy when the 
system is constraint to any OD cavity. Lost cases arise 
in certain OD cavities (for instance, for hard disks, a 
circular cavity of radius larger than R and smaller than 
2i?/v / 3) and thus the FM functional does not recover 
the exact limit for them. This failure of the theory is 
associated to an incorrect low density limit of the third 
order direct correlation function J5| ; in other words, the 
density expansion of the FM functional for hard disks 
or spheres is incorrect already at third order. 

The logic of the construction of FM functional for d- 
dimensional hard spheres requires that the excess free 
energy density be a sum of d terms B3711 

<f d-HS = *£>HS + ^-HS + ' ' ' + *£W (40) 

Further terms constructed on the same logic are identi- 
cally zero. Each of these terms is incorporated starting 
from the first one and trying to compensate for the spu- 
rious terms that arise when two, three, etc., particles are 
incorporated to a OD cavity. When there are lost cases, 
the last one vanishes and thus cannot bring about its 
compensation. The logic of this construction strongly 
relies on the fact that weights have the same shape of 
the particles, as in Rosenfeld's original FM theory OJJ] , 
of which this new constructive method is just a gener- 
alization. 

Applying the same logic to the fluid of aligned hard 
hexagons would lead to the same result and lost cases 
would arise. Yet, there is another method to obtain the 
FM functional for such a fluid, which is the projection 
we have carried out in this work of the fluid of parallel 
hard cubes on a specific plane. The latter fluid does not 
have lost cases because of the particular form of its par- 
ticles (if there is pairwise overlap between three cubes, 
there is necessarily triple overlap as well), and this nice 
property is inherited by the functional for the hexagons. 



As a matter of fact, the resulting functional has a differ- 
ent structure is several respects. First of all, there is a 
weighted density for every one of the maximal OD cav- 
ities conceivable for this system. These include two tri- 
angles and a continuum of irregular hexagons. Because 
of this, weighted densities depend on a parameter that 
gauges the shape of the cavity. This feature is not new: 
it was revealed in the study of FM functionals for lattice 
models ll26i I27I l28i l29i l30ll . But second and more im- 
portantly, there appears a new term, say a "d + 1 term", 
thanks to which the compensation to recover the exact 
OD limit for any OD cavity is guaranteed. As explained 
above, this term cannot be predicted by applying the 
construction logic of the hard disk functional. 

Actually both features are connected: there is an ex- 
tra term because there are several OD maximal cavi- 
ties that contribute weighted densities to the functional. 
What this is telling us is that the functional for hard 
disks or hard spheres is simply incomplete. Cavities 
such as that shown in Fig. [T](and possibly other) should 
make their contribution through new weighted densi- 
ties. Notice in passing that any rotation of the cavity 
of the figure is a new, different, OD cavity and so there 
should also be a continuum of weighted densities, as 
for the hexagons. And accordingly, new terms beyond 
$d_HS are to be expected. How many of them? We 
do not have a definitive answer to this question, but we 
will provide convincing arguments that there will be 
infinitely many. 

The projection we have carried out from three- 
dimensional cubes to aligned hard hexagons can be 
generalized to obtain a FM functional for 2/-gons in 
a straightforward manner. For instance, projections of 
the fourth-dimensional system of hard hypercubes on 
the appropriate plane generates octagons. In general, 
projecting /-dimensional hypercubes on an appropriate 
plane generates aligned hard 2/-gons. Because of the 
structure of the fluid of parallel hard hypercubes lfl2ll 
we know that such functionals will have I terms. Disks 
are obtained as the limit I — ► 00 of 2/-gons, so in this 
limit, the FM functional will be an infinite series. It 
might happen that the series can be summed up and 
produce just a single (more complex) term, but this can 
hardly be foreseen. 

So, should we seek for a FM functional for hard 
disks or hard spheres without lost cases? Well, from the 
arguments above we believe that it would be a pointless 
task, for even if we could overcome the difficulty of its 
construction, we would probably end up with an ex- 
tremely cumbersome functional, useless for practical 
purposes. This does not mean that perhaps including 
some, not all, the missing terms we could obtain im- 
provements on the current functional. This point might 
be worth exploring in the future. 
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APPENDIX A: PROJECTION OF THE FM 
FUNCTIONAL FOR PARALLEL HARD CUBES ON 
THE PLANE x + y + z = 

1. Dimensional crossover 

Using the projection on the XY plane defined by 
r — > Pr, with P given by (\~5l , 

/^phc[/5]= /drSpHcCK^r)}), (Al) 
where 

Pa (Pr) = J dr' u a (P{r - r'))p(Pr'), (A2) 

with the weights defined in (l9ll- (fT2l i. According to 
(fl4] i. in the projected coordinates 



p(Pr) = p(x)5(z). 



(A3) 



This transformation in the density profile allows us to 
integrate in z' , so the weighted densities become 

p a {Pv) = J dx' w Q (P(x - x', z))p(x'), (A4) 

where we are expressing a three-component vector as 
r = (x, z). From the previous equation one can define 
a new set of weighted densities {n a }, which depend 
on a parameter u that corresponds to the z coordinate 
of the above expression, 

n a (x, u) = J dx' O a (x — x', it)p(x'), (A5) 

with a = 0, 1, 2 and Q, a a scalar or vector function 
given by 



f2 2 (x, u) 
fli(x, u) 

Q (x, u) 



w 3 (P(x,u)), 

W2(P(X,«)), 

wi(P(x,«)), 
w (P(x,m)). 



(A6) 
(A7) 
(A8) 
(A9) 



The next step is to obtain the resulting excess free 
energy after the projection. In order to do that, unnec- 
essary degrees of freedom must be eliminated, which 
in this particular case amounts to integrating the z co- 
ordinate in the PHC functional. Consider the first term 
@. Here the integration is immediate, because fio is 
sum of Dirac's deltas 14211. Thus we obtain 



4> 



(i) 

AHH 



1 



in i 



(A10) 



with tiq^ and defined by ( fT6b and ( TTSb , respec- 
tively. 

For the second and third terms in ([3]), direct substi- 
tution of n a leads to 



(2) 
AHH 



du 



(3) 
AHH 



= ; d 



n (u) ■ m(tt) 
l-n 2 (w) 

ni,l(u)rtl,2(M)"l,3(M) 

\l-n 2 {u)? 



(All) 
(A12) 



Here n a (u) is a shorthand for n a (x, u). Notice that 
the integration limits in the above formulas are deter- 
mined by the support of the weighted density 12%, which 
is a product of Heaviside's step functions. It is easy to 
check that this support is |u| < 3L. Hence, accord- 
ing to the expression for ri2(u), maximal OD cavities 
are recovered (see Fig. [5]l when \u\ < L, but when 
L < \u\ < 3i the corresponding cavities are not maxi- 
mal. The appearance of these non-maximal OD cavities 
is a result of our projection procedure, but the func- 
tional cannot explicitly depend on them because any 
information they provided is, by definition, already ac- 
counted for by the maximal OD cavities. In what fol- 
lows we will explain how to get rid of these spurious 
contributions. 



2. Elimination of spurious terms 

We will show here how to eliminate non-maximal 
OD cavities through integration by parts. To this pur- 
pose we introduce the following notation: the weight 
S~2i can be decomposed as follows (in what follows we 
will omit the dependence on (x, u) for simplicity, in- 
troducing only a dependence on u whenever ambiguity 
might arise) 



n x = n{ +) + n\ 



(-) 



(A13) 



where fi^ has the same expression as tt\ in terms 
of functions r y (, but replacing ((x) with S(x =F L), 

respectively. New weighted densities can also be 
introduced associated to these weights. These densities 
have the properties that vanishes when — 3L < 

u < —L, whereas ' does it when L < u < 3L. 
Let us introduce now the differential operator 



D 



d d d d d 
dx du dy du du 



(A14) 
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Acting on n 2 results in 



Dn 2 = -n\ 



(+) , „(-) 



(A15) 



(3) 

As of now, we will only consider the term ^ahh f° r 
—3L < u < —L, which will be denoted 'I'ahh- Ac- 
cording to ( IA15b . in this interval the integrand of 'I'ahh 
can be written as 

<'{n 2 )n[-}n[-}n[^ = 

\ e (^(i)*> 2 ))^ 2) "tW (A16) 

<ren+ 

where is the set of all the permutations a of 
{1,2,3} with positive signature. Using this identity we 
can integrate by parts to obtain 



(3,1) _ n l,2 + n lJ n lJ + rl l,3' >n i,l 



AHH 

\l J_ 3L du *o'(na) £ ( n S(2)"S(3)) («)• 



24[l-r4 } ] 



<renj 



(A17) 

Notice that in the first term, the components of n[ ' 
without explicit dependence on u refer to those of 
n[ '(x, —L), as in ( [T9| >. The second term in the ex- 
pression above can be simplified taking into account 
the new identity 

Ar(2)»h7*(3) + ^(3)^(2) = 2n 0,a(l)- (A18) 



Recall that no is defined through the weight (IA81 . Be- 
sides, 



£ 



A>-(i) ( n i ± ) (2) n i^- ) (3) 



£ (^(2) n S(3) + £> <r(3)"i ± J(2)) ' 



(A19) 

so substitution of JA18t and (|A19t into (|A17b allows 
us to get 



(3,1) l n l,2 "j" ^1^4^ j" n l,3^ n l,l^ 



AHH 



24[1-t4 } ] 



(A20) 



Moreover, as D-n =0 and $ '(n, 2 )nj ) = D$ (n 2 ), 
a new integration by parts in the integral above yields 



(3,1) __ n l,l n lJ j~ n l,2^ n l,3^ n l,3^ n i,l^ 



AHH 



24[l-4 >] 



7^0 } M 1 - "2 



(A21) 



By symmetry, an identical formula can be obtained 
when L < u < 3L, but with the weighted densities 

ria instead of n a . Doing exactly the same partial 
integration in ( IA1 U for any L < \u\ < 3L and gather- 
ing all the contributions, we arrive at the result 

$AHH = $A°HH - ~ £ "0 ±} Ml - 4 ±} ) + 



(±) (±) . (±) (±) . (±) (±) 



24^ 



1 — n 



(±) 



(A22) 

where the residual contribution $ahh ^ s gi ven by 



<T> 1CS 

^AHH 



1 



8 



n (u) ■ m(u) 



du ■ 

l l-m(u) 



du 



—L 



[l-n 2 (u)] 2 



(A23) 



Because of the range — L < u < L of these integrals, 
Eq. dA22l ) only contains contributions from maximal 
0D cavities, as required. 

A third integration by parts in the residual term will 
allow us to write the functional in a more compact way. 
According to ( IA15I ). whenever \u\ < Lwe can write 



ni = 2n < ^ ) ±Bn 2 . 



(A24) 



Then the last integral in ( IA23t can be expressed in a 
symmetric form combining 



1 f du $ > 2 ) £ f[ (2ng 3 ± D jn2 ) («). 
± i=i 

(A25) 

We now expand the product above and realize three 
things. First, that the product containing three deriva- 
tives vanishes. Second, that cross terms containing one 
or two derivatives can be integrated with the aid of 
dABT l and 

^(2)^5(3) + D *(3) n {%2) = ±n 0Ml)> ( A26 > 

which holds for all |u| < L. These terms cancel the 
first one of (fA23) and generate a boundary term iden- 
tical to the last one in ( 1A221 I. This operation justifies 
the numerical coefficient appearing in (|26| |. Third, the 
remaining term, which contains products of the three 
components of nj , produces d27l i. This way we ar- 
rive at the final form of the functional (I24l-d27]i. 
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